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O ABSTRACT 



I 



O ' Aims. The goal of this work is to put constraints on the strength and structure of the magnetic field in the cluster of galaxies A2382. We 
investigate the relationship between magnetic field and Faraday rotation effects in the cluster, using numerical simulations as a reference for 
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the observed polarization properties. 

Methods. For this purpose we present Very Large Array observations at 20 cm and 6 cm of two polarized radio sources embedded in A2382, 
and we obtained detailed rotation measure images for both of them. We simulated random three-dimensional magnetic field models with 
different power spectra and thus produced synthetic rotation measure images. By comparing our simulations with the observed polarization 
properties of the radio sources, we can determine the strength and the power spectrum of intra-cluster magnetic field fluctuations that best 
reproduce the observations. 

Results. The data are consistent with a power law magnetic field power spectrum with the Kolmogorov index n = 11/3, while the outer scale 
of the magnetic field fluctuations is of the order of 35 kpc. The average magnetic field strength at the cluster center is about 3 fiG and decreases 
in the external region as the square root of the electron gas density. The average magnetic field strength in the central 1 Mpc 3 is about 1 //G. 



> 

in 

(N 



O 

o 



X 



Key words. Galaxies:cluster:general - Galaxies:cluster:individual: A2382 - Magnetic fields - Polarization - (Cosmology:) large-scale structure 
of Universe 



£j ■ 1. Introduction 1982, Vallee et al. 1986, Kim et al. 1990, Kim et al. 1991, 

Clarke et al. 2001) or individual clusters. In the latter case one 

The intra-cluster medium (ICM) in clusters of galaxies is analyses deta0ed high resolution RM images (e g Perley & 

known to possess a magnetic field, but its origin and proper- Taylor ml Taylor & Perley m ^ Feretti et al 1995 Feretti et 

ties are not well known. The existence of magnetic fields can al 1999> G(m)ni et al 2001j Taylor et al 2QQh Eilek & Qwen 

be demonstrated with different methods of analysis (see e.g. 2QQ2 ^ Pollack et al 2005j Govoni et al 2QQ6) These data ^ 

the review by Govoni & Feretti 2004, Carilli & Taylor 2002, usually consistent with central magne tic field strengths of a few 

and references therein). The strongest evidence for the pres- ^ Q but stronger fieldfj m found in the inner r£gions of relaxed 

ence of cluster magnetic fields comes from radio observations. cooling core clusteis ^ and can reach values of 10 _ 40 ^ (see 

Magnetic fields are studied through the synchrotron emission e g Taylor £t al 2002) Both for interacting and relaxed clus . 

of cluster-wide diffuse sources, and from studies of the Faraday ters the distribution of radio galaxies is gene rally patchy, 

rotation of polarized radio galaxies. The magnetized plasma indicating that cluster magnetic fie i ds have structures on scales 

that is present between an observer and a radio source changes as sma jj as i q kpc or less 

the properties of the polarized emission from the radio source. Qn the basifj of the available m imageSi increas i ng att ention 
Therefore, the magnetic field strength can be determined with ifj giyen in the literature t0 me power spec trum of the intra- 
the help of X-ray observations of the hot gas, through the inves- duster magnetic field fluctuations . Several studies (EnBlin & 
ligation of the Faraday Rotation Measure (RM) of radio sources Vogt 2QQ ^ Murgia et al 2QQ4) haye shown that detailed m 
located inside or behind the cluster. images of radio galaxies can be used to infer not only the clus- 
RM studies of radio galaxies in clusters have been carried t er magnetic field strength, but also the cluster magnetic field 
out using either statistical samples (e.g. Lawler & Dennison powe r spectrum. The analysis of Vogt & EnBlin (2003, 2005) 
suggests that the power spectrum is of the Kolmogorov type 
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Fig. 1. Left: radio contours of galaxies PKS 2149-158 and PKS 2149-158C obtained at 20 cm superimposed on the red DSS2 
image. The radio image has been obtained by combining all the VLA 20 cm data and by averaging the two IFs within the 20 cm 
band. The sensitivity (1 <t\) is 0.052 mJy/beam. The contour levels start at 3<t\ and are scaled by "V2; the restoring FWHM beam 
is 5.3" x 5.3". Right: zoom of the dumbbell system showing the two radio cores and the respective twin jets. The radio contours 
refer to the A array data at 6 cm band. Levels start at 0.06 mJy/beam (3cri) and increase by a/2; the angular resolution of the radio 
image is 0.4" x 0.5". 



if the auto-correlation length of the magnetic field fluctuations 
is of the order of a few kpc. However, Murgia et al. (2004) 
pointed out that shallower magnetic field power spectra are 
possible if the magnetic field fluctuations extend out to sev- 
eral tens of kpc. Recently, Govoni et al. (2006), who used the 
numerical approach developed by Murgia et al. (2004), derived 
the power spectrum of the intra-cluster magnetic field fluctua- 
tions in A2255, and found that the field strength declines from 
the cluster center outwards, with an average field strength of 
about 1.2 yuG over the central 1 Mpc 3 . They also showed that 
in order to explain both the observed RM of the radio galax- 
ies in A2255 and polarization levels of the radio halo present 
in this cluster (Govoni et al. 2005), the maximum scale of the 
magnetic field fluctuations must be of the order of hundreds of 
kpc with a steepening of the power spectrum from the cluster 
center to the periphery. 

In this paper we present Very Large Array (VLAQ) obser- 
vations at 20 cm and 6 cm of the three polarized radio galaxies 
PKS 2149-158 (A and B) and PKS 2149-158C in the cluster 
Abell 2382; the first two (A and B) form a dumb-bell system. 
A2382 is an ideal case to study RM along different lines-of- 
sight because PKS 2149-158 and PKS 2149-158C are extended 
and highly polarized radio galaxies, located respectively at 4.8' 
and 4.3' from the cluster center. Because the radio sources un- 
der investigation are quite extended, both in angular and lin- 
ear size, they are ideal targets for an analysis of the rotation 
measure distribution: detailed RM images can be constructed 
which can serve as the basis of an accurate study of magnetic 
field power spectrum. 



1 The Very Large Array is a facility of the National Science 
Foundation, operated under cooperative agreement by Associated 
Universities, Inc. 



We follow the numerical approach proposed by Murgia et 
al. (2004), i.e. we study the polarization properties of the ra- 
dio galaxies, at the same time making use of the cluster X-ray 
information. 

The paper is organized as follows. In Sect. [2] we discuss 
details about the radio observations and the data reduction. In 
Sect.[3]we present the total intensity and polarization properties 
of the radio galaxies at 20 and 6 cm. We also describe the mor- 
phology of the sources and the discrete features found in the 
total intensity images. In Sect.|4]we present the X-ray environ- 
ment in which the radio galaxies are embedded. In Sect. [5] we 
show the RM images, discuss the results, and discuss the clus- 
ter magnetic field. In Sect. [6] by following the same approach 
as in Murgia et al. (2004), we introduce the multi-scale mag- 
netic field modelling used to determine the intra-cluster mag- 
netic field strength and structure, and we show the results ob- 
tained with different configurations of the magnetic field power 
spectrum slope. Finally, Sect. [7] summarizes our main conclu- 
sions. 

Throughout this paper we assume a ACDM cosmology 
with H = 71 km s -1 Mpc _1 , D, m = 0.3, and Q A = 0.7. At the 
distance of A2382 (z = 0.0618, Struble & Rood 1999), 1 arcsec 
corresponds to 1.17 kpc. 

2. Radio observations and data reduction 

The radio sources PKS 2149-158 (A and B) and PKS 2149- 
158C were observed at the 20 and 6 cm bands, in all VLA 
configurations. The details of the observations are provided in 
Tab. 1 . The observations in the different arrays were made be- 
tween November 1986 and December 1987. 

The flux densities were brought on the scale of Baars et 
al. (1990) using 3C 286 as primary flux density calibrator. The 
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RIGHT ASCENSION (J20O0) 

Fig. 2. Source PKS 2149-158 (A and B) and PKS 2149-158C: total intensity contours and polarization vectors at 1.46 GHz. The 
angular resolution is 5.3" x 5.3". The first contour level is drawn at -3ct[ and the other contour levels start at 3<xi and are spaced 
by a factor of 2. The lines give the orientation of the electric vector position angle (E-field) and are proportional in length to the 
fractional polarization (10" 50%). 



Table 1. Summary of the VLA observations. 



V 


Bandwidth 


Config. 


Duration 


(GHz) 


(MHz) 




(Hours) 


1.45/1.65 


25 


A 


6.6 


1.46/1.66 


25 


B 


6.3 


1.46/1.66 


50 


C 


3.6 


1.46/1.66 


50 


D 


0.9 


4.82/4.87 


25 


A 


6.1 


4.82/4.87 


25 


B 


5.8 


4.83/4.88 


50 


C 


14.6 


4.83/4.88 


50 


D 


11.3 



The pointing position (J2000) is at RA=2l' , 51 m 57" 



andDEC=-15°37'23". 



same calibrator has been used as the absolute reference for the 
electric vector polarization angle. The phase calibrators were 
nearby point sources observed at intervals of about 30 minutes. 

Calibration and imaging were performed with the NRAO 
Astronomical Image Processing System (AIPS), following 
standard procedures. Several cycles of self-calibration and 
CLEAN were applied to remove residual phase variations. 
Instrumental calibration of the polarization leakage terms was 
obtained using the phase calibrators, which were observed over 
a wide range in parallactic angle. The calibration of the abso- 



lute polarization angle was obtained by assuming for the source 
3C 286 a R-L phase difference of 66° at both 20 and 6 cm. 

The (u, v) data at the same frequencies but from different 
configurations were first handled separately and then combined 
to improve uv-coverage and sensitivity. The separate data sets 
were individually flagged, flux calibrated, and polarization cal- 
ibrated before combination. We combined all arrays at both 20 
cm and 6 cm. Each combined data set was then self-calibrated. 

Images of polarized intensity P = (Q 2 + U 2 ) 1 ^ 2 , fractional 
polarization FPOL = P/I and position angle of polarization 
"P = 0.5 tan _1 (i// Q) were derived from the /, Q, and U images. 

3. Total intensity and polarization properties 

Radio images at 20 cm and 6 cm were obtained by combining 
all the VLA arrays and by averaging the two IFs. The left panel 
of Fig.[T]shows the total intensity contours of the 20 cm image, 
which was restored with a FWHM beam of 5"3, overlayed onto 
the optical imaged from the red Palomar Digitized Sky Survey 
2. A zoom over the dumb-bell system is presented on the right 
panel of Fig.Q] where we show the total intensity contours of 
the 6 cm image, which was restored with a FWHM beam of 

a/4 x a/5. 

2 htpp://archive.eso.org/dss/dss 
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Fig. 3. Source PKS 2149-158 (A and B) and PKS 2149-158C: total intensity contours and polarization vectors at 4.88 GHz. The 
angular resolution is 5.3" x 5.3". The first contour level is drawn at -3ct[ and the other contour levels start at 3cri and are spaced 
by a factor of 2. The lines give the orientation of the electric vector position angle (E-field) and are proportional in length to the 
fractional polarization (10" 50%). 



Table 2. Basic properties of PKS 2149-158 and PKS 2149- 
158C. 



PKS 2149-158: 






Position (J2000) radio core A 


21 ''51'" 54'.' 3 


-15°37'28'.'0 


" radio core B 


21''51'"55'.'0 


-15°37'35'.'8 


Redshift 


0.062 




Total flux density at 20 cm 


424 mJy 




Total radio luminosity at 20 cm 


10 24 - 6 W/Hz 




Overall spectral index o^ocm 


0.95 




Radio source largest linear size 


410 kpc 




PKS2149-158C: 






Position (J2000) radio core 


21*51 m 5978 


-15°38'18'.'2 


Redshift 


0.060 




Total flux density at 20 cm 


138 mJy 




Total radio luminosity at 20 cm 


10 240 W/Hz 




Overall spectral index a^oran 


0.86 




Radio source largest linear size 


210 kpc 





We use the convention S v cc y 



The basic properties of PKS2149-158 and PKS2149-158C 
are given in Tab. [2] Flux densities were estimated after having 
applied the primary beam correction, by integrating in the same 
area the surface brightness down to the noise level. 



For the purpose of the polarization and RM analysis, in- 
tensity and polarization images were produced also for each IF 
separately. The relevant information on these images is listed 
in Tab. [3] Total intensity contours and polarization vectors at 
1.46 GHz and 4.88 GHz are shown in Fig s.|2] and [3] respec- 
tively. Vectors represent the orientation of the projected E-field 
and are proportional in length to the fractional polarization. In 
the fractional polarization images FPOL, we included only the 
points with / > 50"!. 

In the following we give a brief description of the individual 
sources. 

3.1. PKS 2149-158 (A and B) 

PKS 2149-158 (FR class I) is a double system composed of two 
nearly equally bright elliptical galaxies in a common envelope 
("dumb-bell" galaxy). The radio source was first mapped by 
Parma et al. (1991) at 1.4 GHz. Both galaxies of the dumb- 
bell system are radio-emitters, forming a double twin jet system 
like 3C75 (Owen et al. 1985). The two radio cores, labelled 
A and B in Fig.Q] are separated in projection by 13'.'5 which 
corresponds to 15.8 kpc. Their position is reported in Tab. [2] 

The radio morphology of PKS 2149-158, which shows reg- 
ular large amplitude oscillations, is rather unusual and can be 
interpreted in terms of two distinct radio sources whose jets 
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Table 3. Total and polarization intensity radio images, for each individual IF, restored with a FWHM beam of 5.3" x 5.3". 



Config. 


V 


Beam 


o"i 


<r Q 


0"U 




(GHz) 


(") 


(mjy/beam) 


(mjy/beam) 


(mjy/beam) 


A+B+C+D 


1.46 


5.3x5.3 


0.054 


0.024 


0.024 


// 


1.66 


// 


0.047 


0.027 


0.027 


// 


4.83 


// 


0.022 


0.019 


0.019 


n 


4.88 


// 


0.020 


0.019 


0.019 



are strongly interacting. The true (three-dimensional) source 
structure is undoubtedly even more complex, because of the ra- 
dio jets may very well bend along the line-of-sight. The max- 
imum angular size of the radio source at 20 cm is about 5.8' 
(=410 kpc). 

The magnetic field configuration, as traced by the 6 cm im- 
ages, is initially transverse in the jets and becomes circumfer- 
ential in the lobes. However, there are systematic differences 
between the linear polarization at 20 cm and 6 cm because of 
the Faraday rotation effect. At the higher frequency (Fig.[3]l, the 
E-vectors in both jets and lobes are rather ordered and the de- 
gree of polarization is high. In contrast, the ^-vectors at 20 cm 
(Fig.|2j in the center of the lobes, have a much more chaotic 
distribution. The mean fractional polarization is =d2% at 20cm 
and ^23% at 6 cm. The signal is affected by beam depolariza- 
tion at longer wavelength. 

3.2. PKS2149-158C 

The radio source (FR class I) is associated with a single ellip- 
tical galaxy and is unrelated to the dumb-bell system. Its total 
extent at 20cm is about 3' (- 210kpc). The source has a dou- 
ble asymmetric structure with the two jets emanating from the 
core to the north-east and south-west. The north jet and south 
jet extend out to 9.2" and 18" respectively from the core. Of 
course, the north-east jet may appear distorted by projection 
effects and probably bends along the line-of-sight. 

Figs. |2] and [3] show the polarization images of the source 
at 20 cm and 6 cm. The magnetic field configuration, as traced 
by the 6 cm images, is initially transverse to the jets, parallel 
to the southern lobe, and circumferential in the northern lobe. 
The mean fractional polarization is 7% at 20 cm and ^25% 
at 6 cm. As for the dumb-bell system, the signal is affected by 
beam depolarization at longer wavelength. 

4. X-ray environment 

The cluster A2382 has been observed in X-rays with the Rosat 
satellite. The left panel of Fig. [4] shows total intensity con- 
tours at 6 cm superposed on the Rosat PSPC archive image 
(800227p) of A2382. The X-ray image represents intensity 
in the 0.1-2.4keV band. It has been corrected for the back- 
ground, divided by the exposure map (a 17ksec exposure) 
and smoothed with a Gaussian of cr — 30". The centroid of 
the image is located at RA=21 /! 51 m 55"DEC=-15°42'26"; the 
X-ray emission extends up to more than 15'. The radio-X over- 
lay shows that the two radio sources PKS 2149-158 and PKS 
2149-158C are offset to the north of the cluster centre by about 
4.8' (340 kpc) and 4.3' (300 kpc) respectively. 



4.1. X-ray surface brightness profile 

The right panel of Fig. [4] shows the X-ray surface brightness 
(Sx) profile of A2382. The profile was obtained by aver- 
aging the 0. 1-2.4 keV Rosat image (corrected for the back- 
ground and divided by the exposure) in concentric annuli of 
30" in size, centered on the X-ray centroid. Point sources 
have been masked. We converted the X-ray surface brightness 
from counts/skypixel/sec to erg cm^sec" 1 sterad" 1 by using the 
PIMMS0 software (Mukai 2007). In this conversion the X-ray 
emission was approximated by a Raymond-Smith model with 
a mean cluster temperature kT^2.9keV (Ebeling et al. 1996), 
metal solar abundance Z-0.6 and photoelectric absorption col- 
umn density nH=4.1xl0 20 ctrr 2 . 

Because a strong emission core is present in the inner 
100 kpc of the cluster, the observed Sx cannot be described by 
the simple y6-model (Cavaliere & Fusco-Fermiano, 1976) for 
the gas density: 

n e (r) = n (l+r 2 /r 2 c y^, (1) 

where r, no and r c are the distance from the cluster X-ray cen- 
troid, the central electron density, and the cluster core radius, 
respectively. 

Therefore, we first fitted the radial X-ray surface brightness 
only for r > 100 kpc using the y6-model with the three free 
parameters /?ext> r c EXT , n o EXT , resulting in a best fit withx EXT 
1 1 .8. for 16 degrees of freedom. We then used these values in 
a subsequent fit of the S x profile by combining a double 
model: 

n e (r) = « ()lNT (l + r 2 /4 T )-iANx + „ fc(1 + f/^yfan (2) 

The double /3-model is based on 6 parameters (VJint, r ClNT , «0int> 
/3ext, '"cext' n o EXT ) °f which the last three are fixed at the values 
calculated by the previous fit. In this case, the final x 2 NT is 15.0 
for 19 degrees of freedom. The fits of the two y6-model (single 
and double) are shown respectively as dashed and solid line in 
Fig.|4](right panel). The best fit parameters are listed in Tab. [4] 
Overall the model fit is very good. It gives a central gas density 
of «o = "Oint + w o EXT - 5 x 10~ 3 crrT 3 and a outer core radius 
of 373 kpc with a outer /3ext of 0.9. The inner yS INT > 0.7 
and r CmT > 1 1 kpc are lower limits. This means that larger and 
larger values of the inner core radius and /3 still satisfy the data, 
provided that these two parameters grow together. 

5. Rotation measure images 

Polarized radiation from cluster and background radio galax- 
ies may be rotated by the Faraday effect if magnetic fields are 

3 http://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html 
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Fig. 4. Left: radio contours of galaxies PKS 2149-158 and PKS 2149-158C at 6 cm superposed on the Rosat PSPC X-ray 
image. The radio image has been obtained by combining all the VLA arrays and by averaging the two IFs of the 6 cm band. 
The sensitivity (1 o~i) is 0.015 mJy/beam. The contour levels start at 3cn and are scaled by V2; the restoring FWHM beam is 
5'.'3 x 5'.'3'. Right: surface brightness profile of A2383 in the band 0. 1-2.4 keV band. The dashed and the solid lines represent the 
best fit of, respectively, the single y6-model and the double y6-model described in the text. 



present in the intra-cluster medium. Linearly polarized elec- 
tromagnetic radiation passing through a magnetized ionized 
medium suffers a rotation of the plane of polarization: 

y 0bs {X) = x ¥int + {X) 2 xRM (3) 

where ^obsiA) is the position angle observed at a wavelength 
A, is the intrinsic position angle and the Rotation Measure 
(RM) is related to the electron density (n e ), the magnetic field 
along the line-of-sight (By), and the path-length (L) through the 
intracluster medium according to: 

n e [cnr 3 ]B\\ [fiG]dl (4) 



The position angle of the plane of polarization is an ob- 
servable quantity, therefore, images of rotation measure can be 
obtained by a linear fit of the polarization angle as a function 
of A 2 (see e.g. AIPS task RM or the PACERMAN algorithm 

Table 4. A2382 X-ray surface brightness profile best fitting pa- 
rameters. 



Parameter 


Value 


l-o" range 


Units 


^°EXT 


1.2 


1.0-1.5 


10 3 cm" 3 


*" C EXT 


373 


244-625 


kpc 


fiEXT 


0.9 


0.7-1.6 




^Ixt/NDF 


0.74(11.8/16) 






"o, NT 


3.8 


1.8- 20 


10 3 cm 3 


y 

'INT 


65.7 


> 11 


kpc 


/?INT 


1.7 


> 0.7 




*? NT /NDF 


0.79(15.0/19) 







by Dolag et al. 2005). As is well known, determination of the 
rotation measure is complicated because of vm ambiguities in 
the observed ^obs- Removal of these ambiguities requires ob- 
servations at many frequencies that are well-spaced in A 1 . 

We implemented an RM fit algorithm in the FARADAY 
tool (Murgia et al. 2004). Given the U and Q maps at each fre- 
quency as inputs the task UQjoJtM produces the RM and the 
intrinsic polarization angle, both with relative error maps, and 
a x 1 map, obtained by fitting the observed polarization angle 
images. In order to improve the RM image, the algorithm can 
be iterated several times, by using the RM information in the 
high signal-to-noise region and thus improving computations 
in adjacent bad pixels. 

Fig. [5] shows the image of the rotation measure of PKS 
2149-158 and PKS 2149-158C computed using the polariza- 
tion Q and U maps at the frequencies 1.46, 1.66, 4.83 and 
4.88 GHz, with a resolution of 5'.'3. Contours refer to the total 
intensity image at 6 cm. The RM was calculated only in pixels 
with I>3cr! at 6 cm. In the RM image thus obtained we blanked 
all pixels with an fitting error greater than 10rad/m 2 . The bulk 
of RM values range from about -100 rad/m 2 up to 100 rad/m 2 . 
Alternating positive and negative RM patches are apparent in 
each source with RM fluctuations down to scales of a few kpc. 

In order to quantify the goodness of the fit of the rotation 
measure image, we show, in Fig. [6] the RM error image and 
the x 1 map. The average error of the RM as given by the fit 
procedure is about 5 rad/m 2 while the average reduced x 1 is 
1.4. 

We can characterize the RM distribution in terms 
of a mean ((RM) ) and root mean square (ctrm )■ 
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Fig. 5. Images of the rotation measure computed using the polarization Q and U maps at the frequencies 1 .46, 1 .66, 4.83 and 4.88 
GHz with a resolution of 5'.'3 x 5'.'3. Contours refer to the total intensity image at 6 cm. The sensitivity (1 <xr) is 0.015 mJy/beam; 
the contour levels start at 3o"i and are scaled by a factor of 2. 



Fig-El shows the histograms of the RM distribution for 
the two sources. The distributions of RM have approxi- 
mately zero mean; in detail: for PKS 2149-158 we found 
(RM) =0.8rad/m 2 , and ctrm =31 rad/m 2 , while for PKS 2149- 
158C (RM) =-4.2 rad/m 2 , and cr RM =46 rad/m 2 . 

In order to verify the polarization angle linearity with A 2 , in 
the RM map we selected some pixels corresponding to source 
regions with high and low RM and^- 2 values. Fig. [8] shows the 
fits computed in such positions, indicated in the inserted image. 
The data are well represented by a linear A 2 relation. These 
results are in agreement with the interpretation that external 
Faraday rotation is the dominant mechanism in the sources, al- 
though a much larger A 2 range, measured with a finer sampling, 
would be needed to confirm this hypothesis unambiguously. 

Tab. [5] reports, for the two sources (both separate and to- 
gether), the (RM) the ct^m and the maximum (|RM max |) ab- 



solute value of the RM distribution. These data were not cor- 
rected for the Galactic contribution, which is probably neg- 
ligible. In fact, in galactic coordinates A2382 is located at 
lon=38.91° and lat=-46.93° and based on the average RM 
for extragalactic sources published by Simard-Normandin et 
al. (1981), the RM Galactic contribution in a region of about 
10° centered on A2382 is expected to be about -5 rad/m 2 . 

The RM results are consistent with the interpretation that 
the external Faraday screen is mostly due to the intracluster 
medium, indeed the source located in projection nearest to the 
cluster centre, i.e. PKS 2149-158C, also has a higher RM. 

The RM structures on small scales can be explained by the 
fact that the cluster magnetic field fluctuates on scales smaller 
than the size of the sources. These results suggest that it is nec- 
essary to consider a cluster magnetic field that fluctuates over a 
wide range of spatial scales. 
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Fig. 6. Left: RM error image. Right: ^-reduced map of the fit computed to obtain RM image. Contours refer to the total intensity 
image at 6 cm. The sensitivity (1 o^) is 0.015 mJy/beam; the contour levels start at 3cri and are scaled by a factor of 2. 




RM (rad m" 1 ) RM (rad m- 8 ) 



Fig. 7. Histograms of the rotation measure images of PKS 2149-158 (left) and PKS 2148-158C (right). 



Table 5. Rotation measure values. 



Source 


Distance 


< RM > 


C RM 


|RM max | 




(kpc) 


(rad/m 2 ) 


(rad/m 2 ) 


(rad/m 2 ) 


PKS 2149-158 


340 


0.8 


31 


150 


PKS2149-158C 


300 


-4.2 


46 


177 


Both sources 




-0.5 


35 


177 



Col. 1: Source; Col. 2: Projected distance from the X-ray centroid; 
Col. 3: Mean of the RM distribution; Col. 4: RMS of the 
RM distribution; Col. 5: Maximum absolute value of the 
RM distribution. 



6. RM analysis: characterization of the intracluster 
magnetic field power spectrum 

The software package FARADAY (Murgia et al. 2004) permits 
the study of cluster magnetic fields by comparing the observed 
RM with simulated RM images obtained by considering three- 
dimensional multi-scale cluster magnetic field models. In fact, 
given a three-dimensional magnetic field model and the den- 
sity distribution of the intra-cluster gas, FARADAY calculates 
the expected RM image by integrating Eq.|4]numerically. In the 
specific case of A2382, the integration is performed from the 
cluster center up to three core radii (~ 1.1 Mpc) along the line- 
of-sight, i.e. both sources are supposed to lie in a plane which is 
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A s (m') A= (m 2 ) A E (m*J 




A* (m*) A" (m 2 ) A' (m*) 



Fig. 8. Sample plots of the E-vector position angle ^obs against A 2 at different locations in PKS 2149-158 (top panels) and PKS 
2148-158C (bottom panels). The exact position of the points in the sources is shown in the insets of the individual panels. The 
solid line represent the best fit of the A 2 -law to the data. 



perpendicular to line-of-sight and intercepts the cluster centre. 
In the following we neglect the three-dimensional structure of 
the radio sources and assume that all the Faraday rotation oc- 
curs in the intracluster medium in between us and the sources. 



The software has been adapted in such a way that it is now 
possible to treat the simulations in the same way as the obser- 
vations, i.e., using the same fit procedure that is used to derive 
the observed RM. In particular, by using a source model of the 
intrinsic fractional polarization, polarization angle, and the ob- 
served I, U, and Q images, the task RMJoJJQ produces, at 
each frequency and with the same noise as the data, the ex- 
pected Q and U images corresponding to the simulated RM. 
Furthermore, we can take into account both the beam and band- 
width depolarization effects. The synthetic U and Q images can 
be then processed by the task UQjoJtM, resulting in final sim- 
ulated RM images which are filtered with the same algorithm 
as the observations. In this way, the properties of the noise in 
the simulated RM images are very close to those of the data. 
In particular, non-linear effects produced by the fit procedure, 
such as those due to the «7r-ambiguities, are included in the 
simulations as well. 



6.1. The magnetic field model 

We considered the power spectrum of the cluster magnetic field 
to be a power law of the typeQ 

\B k \ 2 (x r n (5) 

in a three-dimensional cubical box. The simulations begin in 
the Fourier space by extracting the amplitude of the magnetic 
field fluctuations from a Rayleigh distribution whose standard 
deviation varies with the wave number according to Eq.0 The 
phase of the magnetic field fluctuations is random. We then per- 
formed a three-dimensional Fast Fourier Transform (FFT) in- 
version to produce the magnetic field in the real space domain. 
The field is Gaussian and isotropic, in the sense that there is no 
privileged direction in space for the magnetic field fluctuations. 
We used a grid size of 1024 3 pixels in size in the wave number 
(Fourier) domain, which allows us to study the magnetic field 
fluctuations over a range of spatial scales more than two-order 
of magnitudes wideH 

4 Note that throughout this paper the power spectra are expressed 
as vectorial forms in fc-space. The one-dimensional forms can be 
obtained by multiplying by 4nk 2 and 2nk respectively the three 
and two-dimensional power spectra. According to this notation the 
Kolmogorov spectral index is n = 11/3. 

5 Here we refer to the length A as the magnetic field reversal scale. 
In this way, A corresponds to a half-wavelength, i.e. A = 0.5 ■ (2n/k). 
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Table 6. Magnetic field model parameters. 

B Average magnetic field strength at the cluster centre 

n Power spectrum spectral index; |6 t | 2 oc ]c" 

Amin Minimum scale of the magnetic field fluctuations 

A max Outer scale of the magnetic field fluctuations 

rj Magnetic field radial profile slope; <B)(r) = (Bp) [^j 



In this set of simulations, the slope and the range of spa- 
tial scales of the magnetic field fluctuation is the same all over 
the cluster volume. However, the normalization of the power 
spectrum decreases with the distance from the cluster center. 
In particular, the average magnetic field strength varies accord- 
ing to: 



<B>(r) = (B Q > 



n e (r) 



no 



(6) 



where (Bq) is the average magnetic field strength at the 
cluster centre while n e (r) is the thermal electron gas density 
assumed to follow the double /3-model profile described in 
Sect.O 

The adopted magnetic field model has five free parameters 
(see Tab.©: B , n, A„„„, A max , and 77. 

By varying all these parameters we obtain synthetic RM 
images characterized by very different statistics and structure. 
Our purpose is to find the combination of model parameters 
that gives the best representation of the observed distribution 
of o~ rm and {RM) across the sources and as well as their RM 
auto-correlation function. 

Ideally, one would like to fit all the five free parameters 
simultaneously. However, in our case this is not very practi- 
cal, because of the computational burden caused by the FFT 
inversion. Therefore we performed a series of simulations that 
search the best magnetic field power spectrum by varying at 
most one or two parameters at a time, while keeping the others 
fixed. We found that there are two main degeneracies between 
the model parameters. The first one is between n and A max : the 
higher is n the lower is A max . The second one is between 77 and 
Bq: the higher is rj the higher is Bq. This means that different 
combinations of these parameters may yield an equally good fit 
to the data. 

In Sect. l6.2l we show the results obtained first by fixing A ml „ 
and A max while varying n. Then we give the results obtained 
by fixing the spectral index at the Kolmogorov value (n=l 1/3) 
while varying A max . In both cases we considered A„„„ = 6kpc 
and 77=0.5. The choice for 77 is motivated in Sect. 16.31 where 
we also analyze how the 77 parameter affects the magnetic field 
strength. The choice for A m ,„ is supported by observations that 
reveal RM fluctuations on small scales. However, in Sect. 16.41 
we analyze the effect of the magnetic field minimum scale on 
the polarization properties of the observed radio galaxies. 

6.2. RM statistics and auto-correlation function 

In the following we compare the simulated and observed RM 
images. In order to assess if a given magnetic field power 



spectrum is able to reproduce the data, we considered two ap- 
proaches: i) we analyzed the RM statistics (<x RM and (RM)) 
calculated over areas of increasing size, and ii) we compared 
the RM auto-correlation functions. 

In order to calculate the RM statistics, we covered the RM 
images with a regular grid of boxes of a given size. We then 
calculated a global average of all the cr RM and (RM) values 
found in each box. By changing the size of the boxes in the 
grid we obtained a trend of the average <r RM and (RM) as a 
function of the box size. We varied the size of the boxes from 
a minimum size of 15 kpc (49 boxes) up to a maximum size of 
300kpc(l box). 

The RM auto-correlation function is calculated as 



A(r) =< RM(x,y) ■ RM(x + dx,y + dy) > 



(7) 



where r = yjdx 2 + dy 2 , while the average is taken over all the 
positions (x,y) in the RM images, excluding blanked pixels. It 
is worth mentioning that A(0) = (RM 2 ) - cr 2 RM + (RM) 2 . 

In Fig. [9] we show the results of a set of simulations ob- 
tained by fixing 77 = 0.5, A m ,„ = 6 kpc, A max = 128 kpc, and by 
varying Bo for three different values of n = 1,2, 3. The choice 
of 77 and Amm is discussed in Sections I6.3l and l6~4l respectively. 
The value of A max has been arbitrarily fixed at 128 kpc. This 
choice is motivated by the evidence of a zero (RM) in both ra- 
dio galaxies, which indicates that the largest scales of the mag- 
netic field fluctuations are smaller than the source size. In the 
top and middle panels of Fig.|9]we present the observed and the 
simulated RM images, using the same colour scale, cellsize and 
resolution. In the bottom left and right panels we show the RM 
statistics and the RM auto-correlation functions, respectively. 

The global ctrm calculated over sources is our most reli- 
able statistical indicator, since it is based on a large number of 
independent measurements. Thus, in the fit procedure we first 
attempt to reproduce the cr RM of the largest box in the statis- 
tics (the 300 kpc dot) by adjusting (Bq) for the three power 
spectra. We obtained a good fit of the cr RM trend for all values 
of n with a central magnetic field strength (Bq) in the range 
3.2 - 4.7 yuG. However, it is clear that n — 3, where most of the 
power of the RM fluctuations are concentrated on A max , results 
in a (RM) level that is much higher than observed, while n = 1, 
where the strongest RM fluctuations are on A mm , generates an 
(RM) lower than observed. The analysis of the (RM) trend sug- 
gests that if A max = 128 kpc the power spectrum spectral index 
should be close to n — 2. This is quantitatively confirmed by 
the values of the reduced x 1 reported in bottom left panel of 
Fig-El Following Govoni et al. (2006), the x 2 has been calcu- 
lated according to 



X' 



-z 



(|(RM obs )| - |(RM sira )|) 2 



eIT KRM ubs >| + eIT KRM sir 



(8) 



>l 



where the errors in the denominator take into account the 
statistical uncertainties in the data as well as in the simula- 
tions^. 



6 For graphical reasons, the simulated trends in left bottom panels 
of Figs.|9]and [I0]are represented as lines. However, the RM statistics 
of the simulated images, which has been calculated with the same set 
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The same behaviour is seen in the RM auto-correlation 
functions shown in the bottom right panel of Fig. [9] A magnetic 
field power spectrum characterized by a spectral index n — 3 
and A max = 128 kpc has too much power at all scales, com- 
pared to the data. On the other hand, a magnetic field power 
spectrum characterized by a spectral index n — 1 generates an 
RM image whose auto-correlation function lies below the ob- 
served one over most of the considered range of scales. The 
RM auto-correlation function corresponding to the intermedi- 
ate case n — 2 gives a better description of the data, confirming 
the result found on the RM statistics. 

Even if the data can be quite well explained by a fiat (n = 2) 
and a broad (A=6 - 128 kpc) magnetic field power spectrum, 
because of the degeneracy existing between n and A max the ob- 
served RM can be explained also by a narrower and steeper 
power spectrum. We produced a second set of simulations in 
which we kept the slope of the power spectrum fixed at the 
Kolmogorov value, n = 11/3, and let (Bo) vary, for three dif- 
ferent values of A max = 25,35,50 kpc. The values of 77 and 
A m i„ are the same as in the previous set of simulations. 

In Fig.[TOJwe compare the observed and simulated RM im- 
ages corresponding to n = 11/3. We found that (Bo) falls in 
the range 2.7 - 4.6 pG, a result which is very close to what 
was found previously. However in this case the degree of sim- 
ilarity between the simulated and the observed RM images 
remarkable. As can be seen in the bottom panel of Fig.fTOl 
the Kolmogorov models reproduce the data better than the 
wider and flatter power spectra considered above. In particu- 
lar, the magnetic field power spectrum characterized by n — 
11/3 and A„ MA = 35 kpc provides an excellent fit to the ob- 
served (RM) profile, yielding a reduced x 2 close to unity. The 
Kolmogorov model with A max = 25 kpc does not have enough 
power on large scales to reproduce the observed (RM) levels. 
The model with n = 1 1/3 and A max = 50 kpc provides a good 
fit to the (RM) statistics but fails to reproduce the observed 
ctrm values on scales below 20 kpc. 

This is further confirmed by the analysis of the RM auto- 
correlation functions shown in the right panel of Fig. [TO] The 
RM auto-correlation function of the Kolmogorov power spec- 
trum with A mnx = 35 kpc is very similar to the observed 
one. The Kolmogorov power spectrum with A„ MA = 25 kpc 
has too much power on small scales and its auto-correlation 
function cuts off faster than the observed one. The model with 
A max = 50 kpc has too much power below 20 kpc and cuts off 
too late in terms of spatial scales, thus failing to reproduce the 
observed RM auto-correlation function. 

To summarize, the analysis of the RM statistics and auto- 
correlation functions reveals that the best fit to the data is ob- 
tained by a Kolmogorov power spectrum with A max = 35 kpc 
and (Bq) = 3.3 pG. 



of boxes used for the data, is affected by the same uncertainty. Two 
effects contribute in determining the amplitude of the scatter. The first 
one is due to the error measurements while the second one is due to the 
statistical variance due to the random nature of the magnetic field in 
the Faraday screen. In our case, the latter is dominant over the former. 



6.3. The magnetic field strength radial profile 

The amount of RM depends on the integral of the product of the 
field intensity and the electron density along the line-of-sight 
(see Eq.@). This dependency results in a degeneracy between 
the 77 and (Bo) parameters. Although 77 does not dramatically 
affect the estimate of the magnetic field power spectrum spec- 
tral index, it can however strongly affect the estimate of the 
magnetic field strength. In particular, the steeper the magnetic 
field radial trend, i.e. the higher is 77, the higher should be (Bo) 
in order to reproduce the observed RM levels. 

Here we justify the choice 77 = 0.5 adopted in Sect 16.21 and 
we critically discuss the estimate of the magnetic field central 
strength as a function of the 77 parameter. 

It is possible to obtain a constraint on the index 77 and there- 
fore allowing measurement of the magnetic field distribution in 
A2382, by comparing the observed and simulated ctrm profile 
as function of distance from the cluster centre. 

In Fig.QT]we show the observed cr RM profile compared with 
the simulations obtained for different values of the 77 parameter, 
for the model which best reproduces the observed RM image 
statistics and auto-correlation function, i.e. n = 11/3, A m ,„ = 6 
kpc, and A nwx = 35 kpc. The radial profiles shown in Fig.fTTl 
have been traced by calculating <x RM in 50 kpc wide concentric 
annuli centred on the cluster X-ray centroid. The simulated ra- 
dial profiles were extended to both smaller and larger distances 
from the centre compared to the range of distances covered by 
PKS 2149-158 and PKS 2148-158C. For this reason, and dif- 
ferent from the procedure described in Sect. |6.2| the simulated 
RM images were not filtered the same way as the observed im- 
ages. Instead, we added in quadrature to the simulated <x RM a 
constant RM noise of 15 rad m~ 2 . This value is the average 
noise introduced by the fit procedure in the case of the consid- 
ered magnetic field model. 

The analysis of the radially average <x RM profiles presented 
in Fig. QT| shows that the best fit of the data is obtained with 
77=0.5 and (Bq) = 3.6 pG. This value of 77 corresponds to a 
magnetic field whose energy density decreases from the clus- 
ter center as the square root of the gas electron density. The 
value of the central magnetic field strength obtained through 
the fit of the radially averaged <x RM profile is indeed in very 
good agreement with that obtained by the fit of the RM statis- 
tics, (Bo) = 3.3 pG. A constant magnetic field, represented by 
the 77 = profile, with a strength of (Bo) = 0.9 pG also provides 
a good fit of the data, although the x 1 is slightly worse than the 
case corresponding to 77=0.5. The steepest magnetic field ra- 
dial profile considered here, 77 = 1, provides the highest central 
magnetic field strength (Bo) = 13 pG but also the worst fit of 
the data. The formal uncertainty in the central magnetic field as 
provided by the fit procedure, i.e. excluding all the systematic 
effects and by keeping fixed all the other four parameters, is of 
the order of (Bo) = 3.6 + 0.5 pG. However, Fig. PTTI shows that 
if we allow both 77 and (Bo) to vary the 1 — cr confidence region 
around the best fit parameters is certainly larger and a safer es- 
timate is to consider the central magnetic field strength in the 
range 1 < (Bo) < 13 pG. 

It should be noted that the above conclusions are based on 
the assumption that PKS 2149-158 and PKS 2148-158C lie at 
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the same distance along the line-of-sight. On the other hand, 
even if different values of the 77 parameter lead to quite differ- 
ent values for the central magnetic field strength, the average 
magnetic field strength over a larger cluster volume is nearly 
the same for the three indices 77. In fact, the average magnetic 
field strength over the central 1 Mpc 3 is almost the same for all 
the three models: (-B)iMpc 3 a 1/jG. 

6.4. The depolarization constraint on the magnetic 
field minimum scale 

In all the models considered so far we fixed the value of the 
minimum scale of magnetic field fluctuations to A m ,„ = 6 kpc. 
This choice was motivated by the fact that this scale matches 
the linear resolution of the radio images, i.e. 6.2 x 6.2 kpc. Here 
we evaluate the effects on the simulated RM images by con- 
sidering a A„„„ that is smaller than the beam. The size of the 
computational box allow us to push the value of A m ,„ down 
to 1 kpc. The most direct consequence arising from a magnetic 
field power spectrum which fluctuates on scale smaller than the 
linear resolution of the observations is the beam depolarization 
effects. With the same FARADAY tool used to filter the sim- 
ulations as well as the observations, it is possible to translate 
an RM image into U and Q images at any given frequency and 
resolution. This allows us to estimate the amount of beam depo- 
larization resulting from different value of A ml „. In Fig.fTSI we 
show the comparison of the observed and simulated fractional 
polarization in the case of the best fit magnetic field model with 
n = 11/3, 77 = 0.5, and A max = 35 kpc. Middle and bottom 
panels of Fig. Q~2] show the simulated fractional polarization for 
A m! -„ = 6 and 1 kpc, respectively. The amount of depolariza- 
tion between 4.88 and 1.46 GHz is not only in good agreement 
with the data, but also almost the same for the two scales. In 
the case of the Kolmogorov spectrum the beam depolarization 
is therefore very weak. This can be explained by the fact that, 
for such a steep spectrum model, most of the magnetic energy 
density resides in the large scales fluctuations. Thus, the sim- 
ulated RM images are almost insensitive to changes on A m! „. 
However, for the same reason we cannot put a lower limit on it. 
It should be considered that lowering A„„„ results in an higher 
magnetic field strength. In the case of the Kolmogorov spec- 
trum with n — 11/3, 77 = 0.5, and A max = 35 kpc, by lowering 
A m ;„ from 6 to 1 kpc requires an increase in the magnetic field 
strength from (Z?o) = 3 to 5 /jG on order to explain the observed 
RM values. This is due to the fact that the magnetic field auto- 
correlation length, Afl , is also smaller and since the RM scales 
as (RM 2 ) oc (B ) ■ A l J 2 (see Eq. 15 in Murgia et al. 2004) we 
need to increase (Bo) if Ag. is lowered. 

The same considerations apply to the model with n — 3 and 
n = 2. The situation is different for n = 1. In this case most of 
the power of the RM fluctuation is concentrated in the small 
scales. Therefore, lowering A m ,„ leads to a significant beam 
depolarization. This is illustrated in Fig. [13] where we report 
the simulated fractional polarization in the case of the shallow 
power spectrum with n = 1, 77 = 0.5, and A mnx = 128 kpc. The 
expected fractional polarization at 1.46 GHz decreases from 
about 11% down to 8% when lowering A„„„ from 6 to 1 kpc. 



Thus for a flat power spectrum with n < 2, A m ,„ should not be 
smaller than the linear resolution of the observations in order 
to prevent the depolarization of the signal at low frequencies. 

7. Summary and conclusions 

In this work we studied the strength and structure of the mag- 
netic field in the cluster of galaxies A2382. Following a numer- 
ical approach we investigated the relationship between mag- 
netic field and Faraday rotation effects in this cluster. For this 
purpose, we presented Very Large Array observations at 20 cm 
and 6 cm of two polarized radio sources embedded in A2382, 
and we obtained detailed rotation measure images for both of 
them. We analyzed the X-ray emission of A2382 observed by 
the ROSAT satellite and derived the radial profile of the elec- 
tron gas density. The observed X-ray surface brightness profile 
cannot be described by a simple /3-model due to the presence 
of a core of strong emission in the inner 100 kpc of the clus- 
ter. A double /3-model provides a better fit of the X-ray surface 
brightness profile. There is indeed the possibility that A2382 is 
a cooling core cluster, although new spectroscopic X-ray ob- 
servations are required in order to confirm the hypothesis. 

We simulated random three-dimensional magnetic field 
models with different power spectra and produced synthetic 
RM images. We filtered the synthetic RM images with the same 
fit procedure used to derive the observed RM images in order 
to ensure a proper comparison of the simulations with the data. 
By comparing our simulations with the observed polarization 
properties of the radio sources, we determined the strength and 
the power spectrum of intra-cluster magnetic field fluctuations 
which best reproduce the observation. 

By assuming that PKS 2149-158 and PKS 2148-158C he 
at the same distance along the line-of-sight and neglecting their 
three-dimensional structure, we conclude that the data are con- 
sistent with a power law magnetic field power spectrum with 
the Kolmogorov index n - 11/3 while the largest scales of the 
magnetic field fluctuations are of the order of 35 kpc. The av- 
erage magnetic field strength at the cluster center is about 3 fiG 
and decreases in the external region as the square root of the 
electron gas density. The average magnetic field strength over 
the central 1 Mpc 3 is about 1 /jG. 
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Fig. 9. Comparison of observed and simulated RM images of the sources PKS 2149-158 and PKS 2148-158C in A2382. The 
observed RM image is shown in the top left panel. The simulated images correspond to the value of A max =128 kpc for the three 
values of spectral index n =1,2, and 3. The model with n = 2 kpc reproduces the observed RM statistics and auto-correlation 
function, which are shown in the bottom panels, best. 
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Fig. 10. Comparison of observed and simulated RM images of the sources PKS 2149-158 and PKS 2148-158C in A2382. The 
observed RM image is shown in the top left panel. The simulated images correspond to the Kolmogorov index n = 11/3 for 
the three values of A max =25, 35, and 50 kpc. The model with A„ MA = 35 kpc reproduces the observed RM statistics and 
auto-correlation function, which are shown in the bottom panels, best. 
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Fig. 11. Left panel: Radially averaged <x SM profile as a function of the projected distance from the cluster center. Solid lines 
represent the trends obtained from the simulation for a magnetic field configuration with n = 11/3 and A„ MA = 35 kpc for 
different values of the rj parameter. Right panel: magnetic field strength profile for the three values of the rj parameter. The 
average magnetic field strength over the central 1 Mpc 3 is about 1 yiG for all the three models. 
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Fig. 12. Observed and simulated depolarization for the Kolmogorov power spectrum. The fractional polarization averages com- 
prise both the sources. 



18 



D. Guidetti et al.: The intracluster magnetic field power spectrum in Abell 2382 



04 QB 



1.46 GHz 










<FPOL>=12% 


Observed 


}] 52 06 71 &2G0 

HIGMTASCEMaC 


315154 31 5143 
ft 43000) 


1.46 GHz 






- 








<FPOL>=11% 




n=1 


A min= 6k P c 


31 53 QQ 315154 31 bl da 
RIGHT A3CH4a.CN 420OQ1 


1.46 GHz 


k 

ft 










-i 




<FPOL>=8% 




n=1 


A min =1 k P C 



4.88 GHz 

it 

; " . 

*< 

<FPOL>=23% A 



Observed 



2 1 52 LM 3153 00 31 51 Si 



4.88 GHz 



■ - - — 1 r 



<FPOL>=19% | 



till 

3152 OQ 315] E 

FIGHT ASCEN3CTH ilTQCO] 



4.88 GHz 



<FPOL>=18% 



n=1 A min =1k P c 



3] £305 31530B 311134 

qiGhTASCeWiCft 4HB0| 



21&30G 315-300 3153 E>4 

H.GHI A=EB*3QN <|2MDJ 



Fig. 13. Observed and simulated depolarization for the shallow power spectrum with n = 1. The fractional polarization averages 
comprise both the sources. 



